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The dynamical structure factor of a Coulomb crystal of 
ions is calculated at arbitrary temperature below the melt- 
ing point taking into account multi-phonon processes in the 
harmonic approximation. In a strongly coupled Coulomb ion 
liquid, the static structure factor is split into two parts, a 
Bragg-diffraction-like one, describing incipient long-range or- 
der structures, and an inelastic part corresponding to thermal 
ion density fluctuations. It is assumed that the diffraction- 
like scattering does not lead to the electron relaxation in the 
liquid phase. This assumption, together with the inclusion 
of multi-phonon processes in the crystalline phase, eliminates 
large discontinuities of the transport coefficients (jumps of the 
thermal and electric conductivities, as well as shear viscosity, 
reported previously) at a melting point. 

PACS numbers: 52.25.Fi, 95.30. Qd, 97.20. Rp, 97.60.Jd 



We consider a strongly coupled Coulomb plasma 
(SCCP) of ions immersed in a nearly uniform charge- 
compensating electron gas. The ions may be disordered 
(liquid phase) or arranged in a crystalline lattice. The 
energetically favorable body-centered cubic (bcc) lattice, 
appears at T > T rn « 172 0, where T = (Zef/{aT) 
is the ion-coupling parameter, T is the temperature, 
a = (47rni/3) -1 / 3 , and rii is the ion number density. 

Many astrophysical objects (interiors of white dwarfs, 
massive stars, and giant planets; envelopes of neutron 
stars) are made of such a plasma. Its kinetic properties 
required for various applications are determined mainly 
by electron-ion (ei) scattering. A general framework for 
calculation of these transport properties has been de- 
scribed in (|]. Numerous calculations (e.g., f^-0]), done 
under additional assumption of strong electron degener- 
acy, predict large (a factor of 3-4) discontinuities of the 
electric and thermal conductivities at the melting point. 
In contrast, the thermodynamic quantities in the liquid 
and solid phases, determined solely by ions, are very sim- 
ilar near T = T m (e.g., |8|j^]). This suggests that prop- 
erties of the ion system serving as a main scatterer for 
electrons should vary smoothly through the melting tran- 
sition. In this Letter, we propose a modification of the 
transport theory which removes large jumps of the trans- 
port coefficients. 

The differential ei scattering rate in a SCCP averaged 
over initial and summed over final electron spin states a 



and a' is 
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where N is the total number of ions, p and p' are the elec- 
tron momenta before and after scattering, respectively, 
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p, Two = 



e is the difference between final 



and initial electron energies, and U^ a i a is the matrix ele- 
ment of the operator of elementary ei interaction. S(q, to) 
is the dynamical structure factor of the plasma, the most 
important quantity of the theory. In the liquid regime, 
p(x, t) is the operator of the charge density in units of 
Z\e\: p(x, t) — ni(x, t) — m, where hi(x, t) is the ion den- 
sity operator and = n e /Z takes account of the com- 
pensating electron background with the electron density 
n e . In the solid regime, p(x, t) = fii(x, t) — <5(x — Ri) 
(where Ri is a lattice vector), i.e. the operator of fluctu- 
ations of the charge density. 

Integrating over x and x' we obtain the structure factor 
of the ion density fluctuations in the solid phase in the 
form 
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where is an ion displacement from R^. Expanding U; 
in the phonon normal coordinates and using the Weyl 
operator identity e A e B = e A+B eJ A ' B ^ 2 , we can decom- 
pose 5 so i(q, t) into the elastic (Bragg) and inelastic parts 
S so \(q, t) — S' sol (q) + S' B ' ol (q,t). The elastic part is easily 
calculated (cf. M): 
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where G is a reciprocal lattice vector, and W = W(q) is 
the Debye- Waller factor, exp(— W) = (exp(iq- Uj)) T , 
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In this case, M is the ion mass, v = (k, s), s — 1,2,3 
enumerates phonon modes, k is a phonon wavevector, 
e„ the polarization vector, lo v the frequency, and n„ = 
(e z " — 1) 1 is the mean number of phonons, z v — Tioj v /T. 
For the lattice types of interest (e.g., bcc), W — r^,q 2 /6, 
where r% — (u 2 )t is the mean-squared ion displacement 
(cf. §,§). 

The Bragg scattering of electrons results in the energy 
band structure of the electron Bloch states, but does not 
contribute to the ei collision integral in the kinetic equa- 
tion Indeed, this scattering occurs at the boundaries 
of the Brillouin zones and translates an electron from one 
zone to another. The transition requires change of the 
electron energy by the value of the interband gap; thus 
another particle must be involved to carry the excess en- 
ergy. 

Therefore only the inelastic part of the structure factor 
contributes to the collision integral. The inelastic part 
can be found by the same technique j^] : 

^2 ^ ~»> \a ijv (n v + 1) + at jv n v \ } , (6) 
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where otij„ = exp[ik • (R^ — R 3 ) — iu) v i\. The summa- 
tion over i,j yields the delta function which removes one 
summation over k (included in the sums over v). Thus 
we have n sums over s and n — 1 sums over k in each nth 
term of Eq. (§). 

Retaining the first term n = 1, we recover the one- 
phonon approximation employed in previous works (e.g., 
[§-§]])■ Our point is that this approximation fails near the 
melting point. In fact, the contribution of the n-phonon 
processes (nth term) at T above the Debye temperature 
can be estimated as (qrT) 2n /n\ ~ (kprT) 2n /nl, where 
kp = (37r 2 n e ) 1 / 3 is the electron Fermi wave number, r\ w 
w_2a 2 /r, u-2 = (tUp/LuDph w 13 is a frequency moment 
for a bcc lattice |Sj], uj v is the ion plasma frequency, and 
(. . .} p h denotes averaging over phonon spectrum (e.g., ||) 
in the harmonic-lattice approximation. For instance, for 
Fe plasma at T w r m we obtain a typical value of n ~ 
(kprp) 2 w 3, which is not small. 

An important difference of astrophysical Coulomb 
crystals from the terrestrial metals is that the umklapp 
processes dominate the scattering, because the equivalent 
radius of the Brillouin zone qs = (Gn 2 ^) 1 ^ 3 is smaller 
than typical momentum transfers q ~ kp. For q > 
qs one can approximately replace X^sll ' e v) 2 f( UJ v) ~ * 
q 2 (/(cLv)) p h p0[ . Then the remaining summations in Eq. 
(^ are done explicitly: 

e 2 ^' ol (q,<) = 
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The static structure factor is defined as S(q) = 
J_°° 5(q, lo) dm = S(q, 0). Thus the inelastic contribu- 
tion to 5(q) in the solid is given by setting t = in 
Eq. (§), which yields S'JJq) = 1 - e~ 2W . 

Let us turn to the liquid phase. Numerical simula- 
tions by different authors show the appearance of incip- 
ient long-range order at T > 1. For example, Schmidt 
et al. [0 observed a shear mode at T > 100 in their 
molecular-dynamics experiment along with the familiar 
longitudinal ion plasmon. We have verified that the spec- 
trum of these modes can be described by the phonon 
spectrum averaged over orientations of a crystal. Al- 
though the long-range order does not persist forever, it 
may be well preserved during typical electron scattering 
time. Thus a temporary electron band structure emerges, 
and an associated elastic scattering does not contribute 
to the conduction (as in solid). This is in line with Ed- 
wards's |l^] argument that one should deal with a local 
disorder "observed" by an electron along its mean free 
path, rather than with the global disorder. Therefore we 
suggest to subtract the elastic contribution from the total 
static structure factor Sn q (q) in the liquid (e.g., p3| , p^ |). 
Since in the liquid an electron couples directly to the ion 
number density, the elastic part must have the form p| 
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Then the "inelastic" part which determines the trans- 
port properties is S{[ q (q) = Sy lq (q) — 5( iq (q). There may 
be various types of periodic structures in this regime, 
but they are very similar and we can use the bcc lat- 
tice. We have checked that the result is almost the same 
for face-centered cubic (fee) and hexagonal close-packed 
(hep) lattices. 

At this stage we need to specify the matrix element 
of the elementary ei interaction Uq. a i a . Assuming the 
Coulomb potential screened by the static polarization of 
ideal, relativistic, strongly degenerate (p w p' ~ hkp) 
electrons, we obtain 
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where ep = [m 2 c 4 +p 2 F c 2 ] 1 / 2 is the electron Fermi energy, 
V is the normalization volume, and e(q) is the electron 
longitudinal static dielectric function Eq |. 

The electric (cr) and thermal (k) conductivities, and 
shear viscosity (n) can be written in the form [0 
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FIG. 1. Dependence of the electron electric (left verti- 
cal scale) and thermal (right vertical scale) conductivities of 
carbon plasma at density 10 4 g cm -3 on temperature (lower 
horizontal scale) or ion coupling parameter (upper horizon- 
tal scale). Dashes are traditional results in ion liquid ||(?]] 
and bcc crystal for T m = 172. Solid line is the present 
multi-phonon calculation in solid, extended artificially to high 
T; dots show the present calculation with reduced structure 
factor in liquid at V < 225. Vertical dotted lines correspond 
to r = 172 and 225. 

are the effective collision frequencies. Here, the effective 
Coulomb logarithms are 
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qo = for the liquid phase and qo = qB for the solid 
phase H, and S a , K (q) arc the effective static structure 
factors. In the liquid regime, we approximate S a ,n(q) by 
S[[ q (q) as described above. In the solid regime, we have 

-i i-+oo r+cc -iut 

S ^ = 2^J duJ J_ dt —e - S "°^^ ^ 



S K ( q ) = S a (q) + (^f-±)*Sj q ). 
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with z = hui/T. The integration over u> can be performed 
analytically. The remaining numerical integration over 
t is then facilitated by shifting the integration path in 
complex plane: t — t' — ih/2T, where t' is real. The final 
result reads: 

1 r +OQ At 

S °(l) = 2j_ -^T-^ 2W K(q,T,t% (17) 
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where x — nt'T/h and 
K(q,T,t') 
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Retaining the term oc q 2 in the expansion of K(q,T, 
we recover the standard one-phonon approximation 

Figures and || show temperature dependence of the 
electric and thermal conductivities for carbon plasma at 
density 10 g cm~ 3 and for iron plasma at 10 8 g cm~ 3 , re- 
spectively, calculated in the Born approximation. In spite 
of large differences in densities and chemical elements, 
the figures are fairly similar. Dashes show the traditional 
results calculated with the full structure factor Sn q (q) in 
ion liquid J^,|?J and in the one-phonon approximation [QJ|] 
for bcc crystal (notice that the results of ref. S for fee 
crystals are in error; actually, they are very similar to 
those for bcc). One can see strong jumps of k and a at 
the melting point T m = 172. Solid lines are the present 
results in the solid phase (including multi-phonon pro- 
cesses), while dots show the present results in the liquid 
obtained using the analytic fits for the static structure 
factor Su q (q) at T < 225 [Q by subtracting the long- 
range correlations (see above). For illustration (as sug- 
gested by H.E. DeWitt), we have extended the improved 
results in liquid by shifting artificially the melting point 
to lower T (to T — 225, considering thus supercooled 
liquid) and the improved results in solid by shifting the 
melting to higher T (lower T, superheated crystal). The 
curves for liquid and solid ions match one another quite 
well, and the jumps at the melting point actually disap- 
pear. In a wide temperature range the improved curves 
for liquid and solid almost coincide. We have verified 
that the same is true for all transport coefficients (in- 
cluding shear viscosity) in wide range of densities for a 
number of chemical elements. Thus one can observe that 
electron transport properties in SCCP of ions appear to 
be fairly insensitive to the state of SCCP (liquid or solid, 
bcc or other crystals) . This should be taken into account 
in various astrophysical applications, for instance, in cal- 
culating the temperature growth from the surface into in- 
terior of the isolated neutron stars important for theories 
of neutron star cooling [plfij] . The thermal conductivity 
of neutron-star envelopes to be reconsidered lies exactly 
in the "sensitivity strip" |l6| which strongly affects the 
temperature profiles and neutron-star cooling. 
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FIG. 2. Same as in Fig. |l| but for iron matter at density 
10 s g cm -3 . 
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